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In this work, exact solutions ol static and spherically symmetric space-times are analyzed in f(R) 
modified theories ol gravity coupled to nonlinear electrodynamics. Firstly, we restrict the metric 
fields to one degree ol Ireedom, considering the specific case ol gtt g rr = — 1. Using the dual P 
formalism ol nonlinear electrodynamics an exact general solution is deduced in terms ol the structural 
lunction Hp. In particular, specific exact solutions to the gravitational field equations are found, 
confirming previous results and new pure electric field solutions are found. Secondly, motivated by 
the existence ol regular electric fields at the center, and allowing lor the case ol gtt grr =£ —1, new 
specific solutions are found. Finally, we outline alternative approaches by considering the specific 
case ol constant curvature, followed by the analysis ol a specific form lor f(R)- 
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I. INTRODUCTION 



A central theme in cosmology is the perplexing fact that the Universe is undergoing an accelerated expansion [l|. 

Several candidates, responsible for this expansion, have been proposed in the literature, in particular, dark energy 

models and modified gravity. Amongst the modified theories of gravity, models generalizing the Einstein-Hilbert 

action have been proposed, where a nonlinear function of the curvature scalar, f(R), is introduced in the action. 

5^ These modified theories of gravity seem to provide a natural gravitational alternative to dark energy, and in addition 

&Q, to allow for a unification of the early-time inflation [2| and late-time cosmic speed-up [3, |j| . These models seem to 

explain the four cosmological phases [5[. They are also very useful in high energy physics, in explaining the hierarchy 

(vq , problem and the unification of GUTs with gravity [g]. The possibility that the galactic dynamics of massive test 

^. ■ particles may be understood without the need for dark matter was also considered in the framework of f(R) gravity 

lO \ models p}- One may also generalize the action by considering an explicit coupling between an arbitrary function of 

the scalar curvature, R, and the Lagrangian density of matter [8|. Note that these couplings imply the violation of 

the equivalence principle |9j, which is highly constrained by solar system tests. 

A fundamental issue extensively addressed in the literature is the viability of the proposed f(R) models [lOl Ull Il3| . 
In this context, it has been argued that most f(R) models proposed so far in the metric formalism violate weak field 
solar system constraints [lJj, although viable models do exis t 111 . 14 L 15 L 16[- The issue of stability [17| also plays an 

o 



important role for the viability of cosmological solutions [T2, LLs HH, [jjj, [2(| . In the context of cosmological structure 



formation observations [21|, it has been argued that the inclusion of inhomogeneities is necessary to distinguish 
between dark energy models and modified theories of gravity, and therefore, the evolution of density perturbations 
and the study of perturbation theory in f(R) gravity is of considerable importance [18. l22l l2a. |24| . 

A great deal of attention has also been paid to the issue of static and spherically symmetric solutions of the 
gravitational field equations in f(R) gravity [25|, [2g, (27|. Solutions in the presence of a perfect fluid were also 
analyzed [281 ]. where it was shown that the pressure and energy density profiles do not uniquely determine f(R). In 
addition to this, it was found that matching the exterior Schwarzschild-de Sitter metric to the interior metric leads 
to additional constraints that severely limit the allowed fluid configurations. An interesting approach in searching 
for exact spherically symmetric solutions in f(R) theories of gravity was explored in |29}, via the Noether Symm etry 
Approach, and a general analytic procedure was developed to deal with the Newtonian limit of f(R) gravity in [301 ]. 
Analytical and numerical solutions of the gravitational field equations for stellar configurations in f(R) gravity theories 
were alsopresented [3l|, [32J, [33[, and the generalized Tolman-Oppenheimer-Volkov equations for these theories were 
derived [311 ]. 

In the context of f(R) modified theories of gravity, it was recently shown that power-law inflation and late-time 
cosmic accelerated expansion can be explained by a modified /(i?)-Maxwell theory [361 ] . due to breaking the conformal 
invariance of the electromagnetic field through a non-minimal gravitational coupling. It is interesting to note that 
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such a coupling may generate large-scale magnetic fields. Motivated by these ideas, we consider in this work f(R) 
gravity coupled to nonlinear electrodynamics, and endeavor to search for exact solutions in a static and spherically 
symmetric set-up. In contrast to a non-minimal gravitational coupling, here conformal invariance is not broken. 

In the context of nonlinear electrodynamics, a specific model was proposed by Born and Infeld in 1934 37] founded 
on a principle of finiteness, namely, that a satisfactory theory should avoid physical quantities to become infinite. 
The Born-Infeld model was inspired mainly to remedy the fact that the standard picture of a point particle possesses 
an infinite self-energy, and consisted on placing an upper limit on the electric field strength and considering a finite 
electron radius. Later, Plebahski explored and presented other examples of nonlinear electrodynamic Lagrangians 
[38|] , and showed that the Born-Infeld theory satisfies physically acceptable requirements. Furthermore, nonlinear 
electrodynamics have recently been revived, mainly because these theories appear as effective theories at different levels 
of string/ M-theory, in particular, in Dp— branes and supersymmetric extensions, and non-Abelian generalizations (see 
Ref. [39j for a review). 

Much interest in nonlinear electrodynamic theories has also been aroused in applications to cosmological models 
[40|, in particular, in explaining the inflationary epoch and the late-time accelerated expansion of the universe [41j . 
It is interesting to note that the first exact regular black hole solution in general relativity was found within nonlinear 
electrodynamics |42| . |43j , where the source is a nonlinear electrodynamic field satisfying the weak energy condition, 
and recovering the Maxwell theory in the weak field limit. In fact, general relativistic static and spherically symmetric 
space-times coupled to nonlinear electrodynamics have been extensively analyzed in the literature: regular magnetic 
black holes and monopoles [44]]; regular electrically charged structures, possessing a regular de Sitter center [4a |: 
traversable wormholes [46( and gravastar solutions [47| . 

Thus, as mentioned above, motivated by recent work on a non-minimal Maxwell-/(i?) gravity model [361 ] . in this 
paper f(R) modified theories of gravity coupled to nonlinear electrodynamics are explored, in the context of static 
and spherically symmetric space-times. This paper is outlined in the following manner: In section [Til the action of 
f(R) gravity coupled to nonlinear electrodynamics is introduced, and the respective gravitational field equations and 
electromagnetic equations are presented. In section IIII[ we restrict the metric fields to one degree of freedom, by 
considering the specific case of g u — —g^r, and using the dual P formalism of nonlinear electrodynamics, we present 
exact solutions in terms of the structural function Hp. Subsequently, in section HVl we investigate the situation where 
the two metric fields are related via a power law in r, introducing additional parameters, and derive new specific 
solutions. In section[Vl we present alternative methods of finding exact solutions, first by considering the specific case 
of constant curvature, then by choosing a form for the /(-R), before we conclude in section fVTl 

II. ACTION AND FIELD EQUATIONS 

Throughout this work, we consider a static and spherically symmetric space-time, in curvature coordinates, given 
by the following line clement 

ds 2 = _ e 2a(r) <ft 2 + e 2/J(r) dr 2 +r 2 (de 2 +sill 2 ed ^2 )i (1) 

where the metric fields a and (3 are both arbitrary functions of r. We use geometrized units, c = G = 1. 
The action describing f(R) gravity coupled to nonlinear electrodynamics is given in the following form 



S = 



Zk 



d 4 x , (2) 



where k — 8n, and f(R) is an arbitrary function of the Ricci scalar R. L(F) is a gauge-invariant electromagnetic 
Lagrangian which depends on a single invariant F given by F = F^F^/4 [38|. As usual the antisymmetric Faraday 
tensor F^ u = A v ^ — A^^ is the electromagnetic field and A^ its potential. In Maxwell theory the Lagrangian takes the 
form L(F) = —F/4ir. Nevertheless, we consider more general choices of electromagnetic Lagrangians. The Lagrangian 
may also be constructed using a second invariant G ~ F ia / f F ,iv ', where the asterisk * denotes the Hodge dual with 
respect to g^y. However, we shall only consider F, as this provides interesting enough results. 

A. Gravitational field equations 

Varying the action with respect to g M „ provides the following gravitational field equation 

IrR^ - -xf ' g^v - VpV,,/^ + g^Ufn = kT^ , (3) 



where f R — df /dR, and the stress-energy tensor of the nonlinear electromagnetic field is given by 

T[j, v = g^v L(F) — F^ a F v a Lp , 



(4) 



with Lp = dL/dF . 

Taking into account the symmetries of the geometry given by the metric (JTJ) , the non-zero compatible terms for the 
electromagnetic field tensor are 



F„„ = 2E(x a ) 5* 5 r J + 2B(x a ) 6% 8f , 



(5) 



such that the only non-zero components are F tr = E(x ti ) and Fqa, = B(x^). Thus, the invariant F takes the following 
form 



F 



e -2(Q+/3) E 2 



B 2 



r 4 sin 2 ( 



Consequently, the stress-energy tensor components are given by 



T\ = T\ = L + e 



rpi) 



T*, 



-2( a +l3) E 2 Li 

B 2 . 



r 4 sin i 



(6) 

(7) 
(8) 



The property T \ = T r r imposes a stringent constraint on the field equations, which will be analyzed further below. 
The contraction of the field equation © yields the trace equation 



f R R-2f + 3a.f R = K T 



(9) 



which shows that the Ricci scalar is a fully dynamical degree of freedom. The trace of the stress-energy tensor, 
T = T%, is given by T = 4(X - FL F ). Note that for the Maxwell limit, with L = -F/(Att) and L(F) = -l/(4?r), 
one readily obtains T = 0, and consequently Eq. ((9J in the Maxwell limit reduces to JrR — 2f + 3 Ofn = 0. 

The trace equation © can be used to simplify the field equations and then keep it as a constraint equation. Thus, 
substituting the trace equation into the field equation © , we end up with the following gravitational field equation 



JrR\ - -ft UrR - Uf R - kT) - V^V./ij 



kT* 



(10) 



Now we can use the properties ((7]) and ([8]) of the electromagnetic stress-energy tensor by subtracting the (rr)-(tt) 
and (90)— (tt) components, which provides the following field equations: 



f'A - (« + 0)'fk-- (« + P)'h 



0. 



(11) 



and 



N(r) 



e -2( a +0) E 2 



B 2 



respectively, where we defined the dimensionless function N(r) as 



N(r) 



r 2 e- 2 ?f R 



2a' 



a' + P'-^ 

JR. 



1-ct 

r 



(12) 



(13) 



The prime stands for the derivative with respect to the radial co-ordinate 7'. It is important to note that Eq. pip 
places a constraint on the metric fields and /#, independently of the form of the electromagnetic Lagrangian. In the 
Einstein limit, /# — 1, Eq. (flTT) leads to (a + (3)' = which we will assume in section ITO1 to explore a specific class of 
solutions. 

Note that with help of Eq. (fTTj) . the following relationship 



□/a = e 



-2/3 



%+(<>/ -(?+*)/'* 



(14) 



2\ e 2/3 -l 

a 1 + - {p - a') - a" + —^ 



and the definition of the curvature scalar, provided from the metric, given by 

R = 2e-^ 
the trace equation © may be expressed as 

/ = f R er^ 



a" + a' (/3' - a 1 ) + -(a' + 5/?') + '—^- + 3 (a 1 + -) & 
r r l \ r J } R 



-T. 



(15) 



(16) 



If a(r) and (3(r) are specified, one can obtain fni r ) from the first gravitational equation (jTTJ) and the curvature scalar 
in a parametric form, R(r), from its definition via the metric. Then, once T is known as a function of r, one may in 
principle obtain f(R) as a function of R from Eq. JT| 



B. Electromagnetic field equations: F representation of nonlinear electrodynamics 

The electromagnetic field equations are given by the following relationships 

(F^L F ).^ = 0, (*FH ;M = 0. 



(17) 



The first equation is obtained by varying the action with respect to the electromagnetic potential A^. The second 
relationship, in turn, is deduced from the Bianchi identities. 

Using the electromagnetic field equation (*F A " / ) = 0, we obtain E = E(r) and B = B{9), and from (F^ L F ) = 0, 
we deduce 



EL i 



q e e 



(«+/3) 



B 



q m sin < 



The electric field E is determined from equations (jl2|) and (|18[> . and is given by 

E(r) 



a+/3 



2nq e 



-N{r)± \ N 2 {r) 



{2Kq e q n 



(18) 



(19) 



Note that independently of N(r) the electric field diverges at the center in the presence of a magnetic field, as in the 
general relativistic case [47j |. Thus, to avoid this problematic feature, in the following analysis we consider either a 
purely electric field or a purely magnetic field. 

The physical fields and the other relevant quantities in the purely electric and the purely magnetic case, respectively, 
are summarized in the following table: 





E(r) 


B{6) 


F(r) 


L F (r) 


purely electric 


Kq e 





1 /iV\ 2 
~2~ \^q~J 


nq 2 e 1 

N r 2 


purely magnetic 





q m sin 6 


2 r 4 


N 2 
«9m 



(20) 



In the purely magnetic case the field equations assume a simpler form, Nr 2 ex Lp, than in the purely electric case, 
where Nr 2 oc 1/L F , and the magnetic fields is independent of the metric fields, contrary to the electric field. Therefore 
the F representation of electrodynamics is more suited for finding purely magnetic solutions which, however, involve 
magnetic monopolcs. 



C. Electromagnetic field equations: Dual P formalism 

As introduced above nonlinear electrodynamics is represented in terms of a nonlinear electrodynamic field, F^ v , 
and its invariants. However, one may introduce a dual representation in terms of an auxiliary field P^ u . This stra tegy 
proved to be extremely useful for deriving exact solutions in general relativity, especially in the electric regime [42|, [43| . 
The dual representation is obtained by the following Legendre transformation 



H = 2FL F ~L. 



(21) 



The structural function if is a functional of the invariant P = P illl P p,v jA,. Then the theory is recast in the P 
representation by the following relations 



P, 



LpF,, 



F, 



fits ~ j-/p-t ' ^ , F iiv — tlpiliv , 

where Hp = dH/dP. The invariant P is given by 

p — _ p p^v — __ 

r - 4 ^>^ - 2 



2PH F 



H. 



LpHp 



(47T)" 



-2(a+/3) p2 _ 



r 4 sin 2 i 



The stress-energy tensor in the dual P formalism is written as 

T„ v = g^{2PH P -H)-P, ia P v a Hp, 

and provides the following non-zero components 

1 



rpT 



H 



r 4 sin 



Pes Hp 



= T% = -H-e-^ a+ ^PlHp. 



(22) 

(23) 

(24) 

(25) 
(26) 



The trace of the stress-energy tensor reads T = —4(H — PHp), so that in the Maxwell limit, H = — P/(4ir) and 
Hp = — 1/(4tt), we have T = 0, which is consistent with the F formalism, as outlined in Section Hi Al 
The electromagnetic field equations now read 



P^ ;/J = 



{Hp*P» v ) = 0. 



(27) 



We emphasize that the tensor F„ v = Hp P^ v is the physically relevant quantity. The P invariant may be deduced 
from Eqs. (|27p in an analogous manner as in the F formalism. In the purely electric case, B = we find 



1e 

2r 4 



(28) 



Due to the fact that it does not depend on the metric fields a and f3, this formalism is attractive to find electric solu- 
tions, as opposed to the usual F representation where purely magnetic solutions are easier to find. The gravitational 
field equation (|12|) now takes the simple form 



r 2 N(r) = -KqlH P {r), 



(29) 



where the function N(r) was defined in Eq. (|13[) and describes the gravity side. Through Eqs. ([20]) in the purely 
electric case we can express the electric field in terms of Hp and P as 



E(r) 



%e a+ ?H P 
r 2 



?a+f3 



V-2PH F 



(30) 



In summary, using the dual P formalism, it is easier to find nonlinear electrodynamic solutions than in the F 
formalism, for the specific case of pure electric fields. We shall consider several specific solutions in the following 
section. 



III. SPECIFIC SOLUTIONS: a(r) = -0(r) 



It is highly non-trivial to find general solutions for the field equations of f(R) modified theories of gravity coupled 
to nonlinear electrodynamics. However, restricting the metric fields to one degree of freedom provides very interesting 
solutions which will be analyzed in this section. In this context, the condition (a + j3)' — imposes a(r) = — /3(r), 
where the constant of integration can safely be absorbed by redefining the time co-ordinate. 

In this specific case Eq. (fir]) implies fp(r) = Ar + B. The Einstein limit is achieved by A — > 0, B — > 1 and so we 
define A := A/ B which represents the departure from Einstein gravity while B can be interpreted as rescaling the 
coupling constants. The second field equation ([25)) now provides the following general solution for the metric field in 
terms of Hp: 



2a(r) 



1 



2r2 



3Ci 



Br z 



Af 



r(f) 



df 



2T(r) 

3r 



9i 

B 



Br z 



Af 



dr 



(31) 



where C\ and C 2 are constants of integration and the function T(r) is defined as 

3 



r(r) - 1 - ~Ar + 3A 2 r 2 - 3A 3 r 3 \u(b[A + l/r]\ . 



The electric field, given in Eq. (|30|) . in this case simply provides 



(32) 



E(r) = %H P (r). 



(33) 



Thus, in principle, by choosing a particular nonlinear electrodynamics theory, by specifying Hp, all the physical 
fields are deduced. Note that in order for the electric field to be finite at the center Hp must be oc r e for small r, 
with e > 2. In the following sections we consider specific choices for Hp and find the respective exact solutions. 



f(R) gravity and Maxwell electrodynamics 



Consider the specific case of f(R) gravity coupled to Maxwell electrodynamics, i.e. H(P) = —P/(4n). The field 
equation ([2^]) provides the following exact solution 



a 2a{r) 



1+AD 



2D 



<il 



3r Br 2 



(l + 2AD)Ar - 2C*i? 



1 



(l + 2AD)\n[B[A+l/r] 



2 ,.2 



A l r 



(34) 



where we defined D — (2Aq 2 — C2V-B, which can be interpreted as an effective mass for the /(i?)-Maxwell case. The 
corresponding electric field is simply 



E(r) 



47T7" 2 



(35) 



and, as expected, diverges at the center. 

Note that the vacuum solution, H(P) = 0, in f(R) gravity, can be immediately obtained by setting q e = in the 
Maxwell solution, Eq. ([34]) . 



3 2a(r) 



C 2 A 2C 2 



B 



3Br 



l .^) Ar . 2Cl ^ 



l-^)':[!J\H I.,] 



2^2 



A l r 



(36) 



An interesting difference to the vacuum solution in general relativity is the term linear in r, and the term with the 
logarithm. Note that the former linear term also arises in the vacuum solutions of conformal Weyl gravity [351 ] . 

In order to obtain the Schwarzschild-de Sitter solution, one sets the following values for the constants: .4 = 0, 
C2 = —3BM and C\ — A/6. This result is similar to the analysis outlined in Ref. [27fl. Note also that the A 7^ 
vacuum solution is not asymptotically flat. An interesting solution is obtained by setting C2 = 0, which yields 



; 2a W = l-Ar-2dr 2 



~ +ln( 13 , A -f J,'/- 



A 2 r 2 



(37) 



This solution has no effective mass term. For positive B it is regular at the center but diverges for large r, independently 
of the constants C\ and A. For negative B it shows the opposite behavior. 



For the specific case of general relativity coupled to Maxwell electrodynamics, i.e. f(R) = R and H(P) = — P/(Att), 
the solution reduces to 



>« = i + ^a + 4-2^ 2 

6r r A 



(38) 



which is simply the Reissner-Nordstrom-de Sitter solution by setting Ci = — 3M and C\ = A/6, as shown above. 
Note that the solution ([55)) is equivalent to considering A — and B = 1 in the solution given by Eq. (f3"4l . 

Clearly it is interesting to try to reconstruct the f(R) theory associated with the solution given in Eq. (fM)) . First 
we calculate the Ricci scalar for the given a(r) which reads in parametric form 



R(r) 



2AD fi4 + RA 2 D 
24Ci + 13.4 2 + 36A 3 D - ^f- + 0Ji + bJilJ + 24 (2d + A 3 D) At 



(1 + Ar) 2 
+6 (4Ci - A 2 ) A 2 r 2 - 12A 2 {1 + 2AD)(1 + Ar) 2 In [B[A + 1/r] 



(39) 



Because of the term oc In (^4 + 1/r), however, this cannot simply be inverted to find r(R). Using the trace equation 
(I16p . with T = for the Maxwell case, we find f(r) in parametric form 



f(r) = 



B 



A(9 + 24C a + 13.4 2 + 36A 3 D) + 6 + 6 - 4 ' + g^!g _ 2^^^^ + 6 (4^ 3 d - ^ 5 ) r 2 



2(1 + Ar) 2 

+AA 2 (1 + 12Ci + 6^ 3 D) r - 12yl 3 (1 + 2^D) (1 + Ar) 2 In (fi [>t + 1/r] 



(40) 



In principle one could find the functional form f(R) from these parametric forms but, as mentioned, R(r) can not 
analytically be inverted to find r(R) and substitute into f(r). 

In figure[T]we plot R(r), f(r), and f(R) for specific values of the constants. B and q e scale the gravitational and 
electromagnetic force while C\ only acts as an overall additive constant, so we set B = 1 = q e and C\ = in the plots. 
The sign of the constant D — (2Aq 2 — C2)/ B influences the sign of R and / close to the center. Thus by looking at 
different values of A with fixed Ci = 2, we see the different behaviors: if D is positive f(R) is not a uniquely defined 
function at large distances r from the center. As a consequence, for f(R) to be well-defined everywhere Ci needs to 
satisfy C 2 > 2.4g 2 . 



R(r) and f(r) for B=l, q=l, C =0, C =2 



f (R) for B=l, q =1, C =0, C =2 



f (r) 


A=0.5 


R(r) 


A=0.5 


f (r) 


A=l 


R(r) 


A=l 


f (r) 


A=1.2 


R(r) 


A=1.2 
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FIG. 1: In the left panel we plot R(r) and f(r), while the right panel shows f(R). We choose the coupling constants to be 
B — 1 = q e and the overall additive constant Ci = 0. The sign of the constant D influences the sign of R and / close to the 
center. For the choice Ci — 2 it reduces to D = 2(A — 1) and by plotting / and R for A £ {0.5, 1, 1.2} we see the different 
behaviors: as soon as D is positive (here A > 1) f(R) is not a uniquely defined function at large distances r from the center. 



2. Generalized Maxwell electrodynamics 

To demonstrate the effect of nonlinear electrodynamics we consider a generalized form of the Maxwell theory 
described by the following structural function 



H(P) = 



P 

-in 



1 



/'- 



1 



■x-^ipf 



(41) 



where fi and 5 are the characteristic parameters of the theory. Note that this choice may physically describe strong 
fields, as the second term is now dominant, i.e. for P ^ 1. This Lagrangian possesses the correct Maxwell limit for 
S > 0, i.e. H ~ — P/(in) for P<1. The relevant quantity Hp is then given by 



H P 



in 



[l + n(-2q* e P) s ] 



(42) 



A particularly interesting and simple example is obtained by setting S — 1/4, such that using Eq. ([28]) . Hp takes 
the form 



H f 



1 

47T 



1 






(43) 



Thus, substituting Eq. (|4"3")l into Eq. (f2T)|) , we finally deduce the following solution: 



e 2Q M = i + AD 

3r 



2V t q 2 e (2-A\qM , 2|< ?e | 3 /x 



2Br 2 



5Pr 3 



'1 + 2_4r>).4r - 2Ci? 



(l + 2^P)ln(P[^+l/r] 



A 2 r 2 . 



where now the effective mass is generalized to 



V 



Aq 2 e (2-A\q e \^-C 2 
B 



(44) 



(45) 



Note that f(R) gravity coupled with Maxwell electromagnetism, i.e. Hp = —i/(4n), follows from the above solution 
in the limit of /i = 0, which simply reduces to Eq. (f3~4")> . Note that we can use the solution (|4"4"1) to write R(r) and 
/(r) in parametric form, and finally, in principle, deduce the functional form f(R). However, as outlined in Section 
IIII 1[ R(r) cannot be analytically inverted to find r(R) and substitute into f(r). In addition to this, we do not write 
out the explicit forms of R(r) and f(r) due to their lengthy character. 

Setting A — and B = 1, which is equivalent to general relativity, Eq. (|44p provides a particularly interesting 
solution given by 



,2a(r) _ 



= 1 



2C 2 
3r 



2|ge| 3 At 

5r 3 



2C ir 2 , 



(46) 



which can also be found from Eq. pip . Note the presence of a term proportional to 1/r 3 , which dominates for low 
values of r. This solution tends to the Maxwell-Einstein limit setting /i = 0. 



IV. NEW SOLUTIONS: a(r) / -j3{r) 



Due to the fact that f(R) gravity has more degrees of freedom compared to Einstein gravity, and also in view of 
Ref. [481 ] . it is very interesting to explore the situation of a ^ — /?. However, without specifying a relation between 
a and (3, a specific nonlinear electrodynamics model, or a specific f(R) theory the equations are not closed and 
therefore analytically intractable. In this section we consider the specific example where the two metric fields satisfy 
the following relationship 



a(r) + fi{r) = ln(fcr^) 



(47) 



where k and t are free parameters. This case is particularly interesting since it allows for regular electric fields at the 
center, as will be shown below. 



From the first field equation (fTTj) we find that fii{r) has the following form 

f R (r) = Ar p+{1) + Br p - {1 \ 
where A and B are constants of integration, and the exponents depend on the parameter £ as 

._(£) = i(l+*±^(*+10) + l) • 



p±( 



(48) 



(49) 



In order for p±{£) to be real it is required that either £ > 2\/6 — 5 or £ < — 2\/6 — 5. In the limit £ = the exponents 
become p+(0) = 1 and p_(0) = such that fn(r) = Ar + B as in the case considered in the previous section. We 
plot p± in figure [5] 





FIG. 2: In the left panel we plot p±, the powers of r in the solution for /ij(r), as function of the parameter £. The right panel 
zooms into the right branch of p±. Note that p~ vanishes for £ = 0, is positive for negative £ and vice versa, while p~ — * —2 
for large £. This behavior makes the solution fn = Br p ~ particularly interesting. 



The electric held now reads 



E(r) = q e kr'- 2 H P {r). 



(50) 



An interesting case is £ = 2 where the electric field is constant in the Maxwell limit, Hp = — l/(47r). 

Let us consider the case A = where fair) is a simple power law which has a well defined Einstein limit for £ = 
and B = 1. For this specific case the second field equation (|2"9"]) can in principle be solved for any general structural 
function Hp(r). We define the constant K — Kq^k 2 /B and solve equation (f2"9"| for £ ^ 1,2, which provides the 
following solution 



.,2n 



k 2 r 2l 



2r 2 



(d+K f fP+ +i - 6 H P (f)di 



P+ - 4 (C 2 +K 



(i-l)(p+-2i-2) p+-4 
while the special cases £ = 1, 2 have to be solved separately. For £ — 1, the solution is 
(61n(r)-\/3-3)fc 2 r 2 2r 2 



2 H P {f)di 



-,2o 



3(v / 3-3) 
and for £ — 2, we find 

1 



d+K J f^- i H P 



(r)dr) 



^- 3 (C 2 +K 



Hp{¥) 



= --k z r 4 + 2r 



3(V3-3) 

(d+K f ln(f)Hp(f)df) - ln(r) (c 2 + K f H P (r)df) 



dr 



(51) 



(52) 



(53) 
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In all cases, the solution is not conformally flat due to the first term. 

An interesting case is f(R) gravity coupled to Maxwell electrodynamics where Hp = — l/(47r). The electric field 
for this case is given by 



E ( r ) = ~^ 



(54) 



where for £ = 0, the classical Coulomb field is recovered. For £ < 2 it diverges at the center. Interestingly the electric 
field is constant for £ — 2, as mentioned before. For i > 2 the electric field vanishes at the center and diverges at 
spatial infinity. 

The metric field is then given in the three cases as: 



1*1,1 
£ = 1 

e = 2 



k 2 r 21 



K r p++?>-3 



2(Cir 2 + C 2 rP+- 2 ) 



. 2 a 



(£ - l)(p + - 2£ - 2) 2tt(£ - 1){ P+ + £ - 5) 
3 - v3 f „ 9 . „ ,/^_ 1 K r 



3(2 - VS) 



d r 2 + C 2 r^- 1 ~ £■ ln(r) + (3 - V3)" 1 

47T - 



p+-4 
Va-i 



r 



K 



_± fc 2 r 4 , « .3 , o„2 



— H + 2r 2 [d-C 2 ln(r)] 

Z7T 



(55) 
(56) 
(57) 



In the first case, the exponents of r are positive in all terms for £ > £ cr ;t, where £ cr n = (5 — VT3)/2 ~ 0.7. See the left 
panel of figure [3] for a comparison of the different exponents. At £ = 1 and £ — 2 the hierarchy of the terms change 
which explains why these are special cases. 

A solution is regular at the origin if the function and all its derivatives are finite at r = 0. We verify that the 
solutions (R)6"l) and (f5"7| are not regular at the origin, although they vanish for r = 0. In case of solution ([B3|) we only 
consider 2<<eN for which the electric field is regular at the origin, cf. Eq. (|M|) . For its derivatives to be finite at 
the origin the following exponents of r 



m -- 


= P+ +£-3 = 


1 - 

2 . 


n 2 - 


= P+-2 = 


1 " 

2 . 



U - 5 + ^/£ (£ + 10) + 1 



£ - 3 + ^£ {£ + 10) + 1 



must be natural numbers, i.e. n\,ri2 G N. 

The metric function g rr must also be regular at the origin. Using a(r) + f3(r) 



(58) 
(59) 



ln(fcr £ ), we have 



with the constants 
C 1 



e 2ti = 



k 2 



Ci + C 2 rP+-'- 3 + C 3 (Cir 2 ^-*) + C 2 rP+- 2 - 2 ') 



(£-l)( P+ -2£-2) 



C 2 



/< 



27r(^-l)(p + +^-5) 



C 3 



P+-4 



(60) 



(61) 



We have considered that I > 2, so that we have to impose C\ — to have regularity at the origin. Furthermore the 
metric function e 2 " and its derivatives need to exist at the origin which means the following exponents of r must be 
natural numbers 



m i = P+ ~ & ~ 3 = 

m 2 = p+-2(£+l) = 



-£-5+^£(£+l0) + l 



-3(£+l) + ^(^+10) + l 



(62) 
(63) 



However, it turns out that mi and m 2 are both negative for all £ G N which contradicts the imposition that mi and 
m 2 be natural numbers. Thus, we conclude that the solution (|55p is not regular at the origin. 

It is evident that for £ = 1 and £ = 2 the metric field goes to negative infinity for large r > rn and thus has to be 
matched to an external vacuum solution at a junction interface at r < rn. This behavior is independent of the signs 
of the constants of integration C\ and C 2 . However, for the solution to be positive for small r, we find C 2 > 0. In the 
case of £ > 2 we find the same behavior of the metric field and the same constraint on C 2 . For £ cl -n < £ < 1 we find 
positive solutions for all r if C 2 < and C\ > 0. If C\ < the solution is negative for small r > 0, if C 2 > it tends 
to negative infinity for large r. See right panel of figure [3) 

We also emphasize that through R(r) and f(r) (not written out explicitly due to their extremely lengthy nature) 
expressed in parametric form, the functional form f(R) may in principle be deduced. However, as outlined in Section 
IIII 11 R(r) cannot be analytically inverted to find r(R) and substitute into f(r). 
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FIG. 3: In the left panel the exponents of r in the four terms of e 2a ^ for (, ^ 1, 2 are plotted. For £ > £ cr i t — 0.7 the exponents 
of r are positive in all terms. Negative £ lead to solutions divergent at the center. At £ = 1 and £ — 2 the hierarchy of the terms 
change which explains why these are special cases. The right panel shows e a ^ r ' for £ = 0.8 for four different combinations of 
signs of the integration constants, while for simplicity and transparency we set K = 1 = k 2 throughout the analysis. 



ALTERNATIVE APPROACHES 



Constant curvature 



An interesting alternative is to consider the specific case of constant curvature R(r) = Rq- Note that in this case fn 
is independent of 7', and for simplicity one may set fii(r) = B. Thus, one verifies that Eq. (jlip yields a(r) = — /3(r), 
so that the curvature scalar is given by 



R 



2e 2a(r) 



[4m' + r 2 a" + 2r V) 2 - e~ 2 " + l] 



For constant curvature, R{r) = Rq, this yields the following solution for a(r) 



2a(r) 



1 



2Ci 2C2 Rq 



r r 2 12 

Substituting the metric field into Eq. (|2"!)j) . one deduces Hp, given by 

C 2 B 



H P {r) = - 



2-Kql 



(64) 



(65) 



(66) 



which reduces to the Maxwell type, i.e. Hp = — l/(47r), by setting the constant of integration C2 = q 2 /(2B). 

For this case, i.e. constant curvature, and taking into account that the Maxwell limit implies T = (see Section 
III Ap . the trace equation ^ imposes the following algebraic relationship 



f R R-2f = 0, 



(67) 



so that the form of f(R) needs to obey this algebraic identity. Thus the metric given by Eq. (|65|) is an exact solution 
for the class of solutions f{R), in the Maxwell limit, that satisfy fp(Ro)Ro — 2f(Ro) = 0. For instance, considering 
the case of f(R) = R — fJ- 4 /R, and using the above trace equation yields Rq = ±\f~ii 2 . The case of f(R) = R + "/ 2 R 2 , 
provides Rg = 0. 
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B. Specific gravity theory: f(R) = R + -y 2 R 2 

Another alternative approach is to consider specific choices for the form of f(R). Consider the specific case of 
f(R) = R + "f 2 R 2 , which for a(r) = —/3(r) implies that Ar + B = 1 + j 2 R, with 7 2 = 2~/ 2 . Substituting the value for 
R(r), provides the following solution: 

Note that this solution is consistent with Eq. (TTTJ) , i.e. f'^ir) — 7 2 i?"(r) = 0. 

Now, substituting this solution in Eq. (fl"2"|) , and finally using the relationship P = —q 2 /(2r 4 ), we reconstruct the 
following nonlinear electrodynamic structural function 



, , C X B A 



2irq 2 e 27^2 



3C 2 \ qe \(-2P)i - 2\q e \i(-2P)i + ^1(-2P)-1 + ^(-2P)" 



(69) 



Note that for A = and C\ — q 2 /(2B) it reduces to the Maxwell type, i.e. H = — P/(4tt). However, for A ^ this 
structural function does not tend to the Maxwell limit for P <§C 1. Therefore it is not a viable nonlinear electrodynamic 
theory. This specific case illustrates the difficulty in finding viable nonlinear electrodynamic theories, i.e. with the 
correct Maxwell limit, by explicitly providing a form for f(R). 

VI. CONCLUSION 

The issue of exact static and spherically symmetric solutions in f(R) modified theories of gravity is an important 
theme, mainly due to the analysis of weak field solar system constraints, and the generalization of exact general 
relativistic solutions to f(R) gravity. In this work we have analyzed exact solutions of static and spherically symmetric 
space-times in f(R) modified theories of gravity coupled to nonlinear electrodynamics. Firstly, the metric fields were 
restricted to one degree of freedom, by considering the specific case of g u — —Qrr ■ Using the dual P formalism of 
nonlinear electrodynamics, an exact general solution was found in terms of the structural function Hp. In particular, 
exact solutions to the gravitational field equations were found, confirming previous results and new pure electric field 
solutions were deduced. Secondly, by allowing two degrees of freedom for the metric fields, and motivated by the 
existence of regular electric fields at the center, new solutions were found. Finally, we have also briefly considered 
alternative approaches by analyzing the specific case of constant curvature and secondly, by considering a specific 
form for f(R). 
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